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ABSTRACT 

We present a comparison between the [Ca,C,N/Fe]-mass relations observed in local 
spheroids and the results of a chemical evolution model which already successfully 
reproduces the [Mg/Fc]-mass and the [Fc/H]-mass relations in these systems. We find 
that the [Ca/Fe]-mass relation is naturally explained by such a model without any 
additional assumption. In particular, the observed under-abundance of Ca with respect 
to Mg can be attributed to the different contributions from supernovae Type la and 
supcrnovae Type II to the nucleosynthesis of these two elements. For C and N, we 
consider new stellar yields that take into account stellar mass loss and rotation. These 
yields have been shown to successfully reproduce the C and N abundances in Milky 
Way metal-poor stars. The use of these new stellar yields produces a good agreement 
between the chemical evolution model predictions and the integrated stellar population 
observations for C. In the case of N, the inclusion of fast rotators and stellar mass-loss 
nucleosynthesis prescriptions improves our predictions for the slope of the [N/Fe] vs. 
a relation, but a zero point discrepancy of 0.3 dex remains. This discrepancy cannot 
be removed, either by increasing the N yields or by assuming a larger amount of fast 
rotators in spheroids, because in both cases this leads to an overproduction of the 
N abundances in the gas phase in these galaxies at high redshift (e.g. the Lyman 
Break Galaxy MS1512 cB-58). This work demonstrates that current stellar yields are 
unable to simultaneously reproduce the large mean stellar [(N/Fc)] ratios inferred from 
integrated spectra of elliptical galaxies in SDSS and the low N abundance measured 
in the gas of high redshift spheroids from absorption lines. However, since chemical 
evolution models for the Milky Way computed with the Geneva stellar yields constitute 
at present the only way to account for the N/O, C/O and 12 C/ 13 C abundance ratios 
observed in very metal-poor halo stars , it seems reasonable to suggest that there may 
be uncertainties in either the inferred stellar or gas-phase N abundances at the level 
of ^0.3 dex. 

Key words: galaxies: elliptical and lenticular, CD - galaxies: abundances - galaxies: 
formation and evolution; galaxies: stellar content - galaxies: individual: MS 1512-cB58 
- stars: rotation 



1 INTRODUCTION 

Abundance ratios can be used to constrain both galaxy for- 
mation scenarios and stellar nucleosynthesis. So far, in stud- 
ies of elliptical galaxies, abundance ratios have been used 
with the former aim. Indeed, increasing evidence has ac- 
cumulated over the past decade that the [Mg/Fe] ratio is 
super-solar in the cores of bright galaxies and increases with 
galactic mass (e.g. Faber et al., 1992, Worthey et al. 1992). 



It is also well established that elliptical galaxies obey the 
mass-metallicity relation, namely the stars of the brighter 
spheroids are more metal rich than those of the less lumi- 
nous ones (Carollo et al., 1993). These two relations, along 
with the evolution of the luminosity function (e.g. Bundy et 
al. 2006) led to the so-called downsizing picture: more mas- 
sive galaxies form faster and earlier than less massive ones 
(Matteucci 1994). 

Unfortunately, due to the fact the we cannot resolve 
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stars in elliptical galaxies and because of the uncertainties 
in modelling line-strength indices as a function of abundance 
ratios, elliptical galaxies have not yet been used to constrain 
stellar nucleosynthesis calculations. Up to now, such studies 
have been restricted to the Milky Way (MW; e.g. Chiappini 
et al. 2003a, b, Cescutti et al. 2006) and its satellites (e.g. 
Lanfranchi et al. 2006, 2008 for dwarf spheroidals; Pompeia 
et al. 2008 for the Large Magellanic Cloud), and Damped 
Lyman Alpha systems (e.g. Dessauges-Zavadsky et al. 2007). 

As an example, metal poor halo stars in the MW are ob- 
served to have high levels of N/O and C/O (Spite et al. 2005, 
Akerman et al. 2004, Fabbian et al. 2009). The observed 
trends in C/O and N/O versus O/H at very low metallici- 
ties cannot be reproduced by chemical evolution models in 
which standard yields are adopted (see Chiappini et al. 2005 
for a detailed discussion). In fact, due to the very low metal- 
licities of such stars, the high levels of N/O would have to 
be the result of nucleosynthesis in metal-poor massive stars, 
suggesting these objects must produce primary nitrogen in 
non-negligible quantified]. It has recently been shown that 
fast stellar rotation is a promising mechanism for producing 
primary nitrogen in metal-poor massive stars (e.g. Meynet 
et al. 2006, Hirschi 2007). Chemical evolution models which 
include the yields of fast rotating models at low metallic- 
ity (Z — 10 -8 by mass) can account for the observations 
in normal metal-poor halo stars (Chiappini et al. 2006a,b, 
2008). 

In principle, there is no reason why such fast rotators 
should occur only in the MW halo. These stars could also 
have left their imprints in DLAs and sub-DLAs (the most 
metal poor DLAs found at present). In fact, Pettini et al. 
(2008) find that these metal-poor high-redshift systems also 
show large C/O and N/O ratios, compatible with what is 
found in the MW halo. Given these successes, it is important 
to determine whether the stellar yields adopted in Chiappini 
et al. (2006a) can also give consistent results for other sys- 
tems whose star formation histories differ from that of the 
solar vicinity and the MW, such as DLAs (Chiappini et al. 
in prep) and elliptical galaxies (this work). 

The large number of galaxies observed in the recent 
Sloan Digital Sky Survey (SDSS; York et al. 2000) make it 
possible to construct very high S/N spectra for typical ellip- 
tical galaxies by stacking together the spectra of many sim- 
ilar objects. The resulting data are of high enough quality 
to enable the study of line-strength ratios for multiple dif- 
ferent elements on a more robust basis (Graves et al. 2007). 
At the same time, Schiavon (2007) has produced a new class 
of stellar population models which take into account the ef- 
fect of various abundance patterns. These models allow us 
to estimate not only the total metallicity, age, and [cv/Fe] of 
galaxies, but also the relative abundances of several other 
elements, such as carbon, nitrogen and calcium. 

Graves et al. (2007) find that many properties of the 
stellar populations in red sequence galaxies (such as age, 
[(Fe/H)] and [(Mg/Fe)]) vary with galaxy velocity dispersion 
(a). In particular, more massive galaxies are typically older, 



1 At least a factor of a few hundred more than what was already 
predicted by the models of Meynet & Macdcr (2002) with ro- 
tational velocities of 300 km s _1 (see Chiappini et al. 2005 for 
details). 



more metal-rich, and more enhanced in Mg, as predicted 
by the chemical evolution models of Matteucci (1994) and 
in agreement with previous work (e.g., Trager et al. 2000b, 
Bernardi et al. 2003, Nelan et al. 2005). Using the newly- 
developed capabilities of the Schiavon (2007) models, they 
further measure [(Ca/Fe)], [(C/Fe)], and [(N/Fe)]. They find 
that Ca scales with Fe and, more surprising, that [(C/Fe)] 
and [(N/Fe)] increase with a more strongly than [(Mg/Fe)] 
does. 

These results are crucial new constraints to study the 
role played by stars of different masses in the chemical en- 
richment of elliptical galaxies. This investigation aims at 
reproducing and understanding these abundance trends by 
means of a self-consistent chemical evolution model, whose 
predictions already match many previous observations. In 
this way we can, for the first time, use several abundance ra- 
tios measured in the stars of elliptical galaxies as constraints 
on nucleosynthesis calculations. 

This paper is organized as follows. We first review the 
observational data and their uncertainties in Section 2, then 
present the main ingredients of our chemical evolution mod- 
els in Section 3. In Section 4, we present our results and then 
summarize our main conclusions in Section 5. 



2 OBSERVATIONS 
2.1 The Data 

Graves et al. (2007) use a sample of ~6,000 galaxies from the 
SDSS Main Galaxy Survey (Strauss et al. 2002) to determine 
the typical stellar population properties of elliptical galaxies 
as a function of a. Their sample compares the stellar popula- 
tions of quiescent elliptical galaxies (those with no emission 
lines in their spectra) to those of elliptical galaxies that host 
emission from low ionization nuclear emission-line regions 
(LINERs). In this work, we are only interested in quiescent 
galaxies. Graves et al. sort galaxies into six a bins, then stack 
the spectra of the galaxies in each bin to obtain very high 
S/N mean spectra. They then measure the strength of mul- 
tiple absorption features in the spectra (using the Lick/IDS 
index system pioneered by Burstein et al. 1984) and use the 
models of Schiavon (2007) to convert the line strengths into 
measured values of mean, luminosity- weighted age, [(Fe/H)], 
[(Mg/Fe)], [(Ca/Fe)], [(C/Fe)], and [(N/Fe)] for the stellar 
populations. 

In agreement with many previous authors, Graves et al. 
find that age, [(Fe/H)], and [(Mg/Fe)] all increase with a 
(see their Table 8 for a detailed comparison with previous 
results). In particular, they find that [(Mg/Fe)] ranges from 
+0.13 to +0.26 as a goes from 100 km s _1 to 250 km s -1 . 
They furthermore find that [(Ca/Fe)] has only very weak 
dependence on a and is nearly solar for all galaxies, with 
values in the range —0.01 < [(Ca/Fe)] < +0.03 for galaxies 
with 100 km s _1 < a < 250 km s" 1 . 

Some of the most interesting results from the Graves 
et al. study, and the ones which motivated this analysis, 
were the observed trends in [(C/Fe)] and [(N/Fe)]. Both of 
these elements show strong enhancement with increasing a 
(see Figures 1 and 2). 

An overabundance of C and N in more massive ellip- 
ticals has been suggested by previous authors, but strong 
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evidence has been lacking and there has been considerable 
confusion in the literature. The first hints came from Trager 
et al. (2000a) 's analysis of a sample of local ellipticals. They 
suggested that, in order to best fit their absorption line mea- 
surements for elliptical galaxies, C should be enhanced along 
with Mg and possibly O. Trager et al. therefore included C 
(as well as N) in the group of "enhanced" elements when 
constructing their stellar population models but did not at- 
tempt to quantitatively measure [(C/Fe)] or [(N/Fe)]. Later, 
Sanchez-Blazquez et al (2003) found a positive correlation 
between the strength of a CN absorption band and a, which 
they interpreted as evidence for further enhancement in C 
and N with increasing a. They did not attempt a quantita- 
tive conversion into abundance ratios. 

Sanchez-Blazquez et al. further presented some evi- 
dence that galaxies in denser environments were somewhat 
less enriched in C and N with respect to their low density 
counterparts. However, Toloba et al. (2009) have recently 
claimed that this correlation disappears if one uses an NH 
feature at shorter wavelength. To further complicate mat- 
ters, Clemens et al. (2006) corroborate the environment- 
correlation proposed by Sanchez-Blazquez et al., but find 
an anti-correlation of [(C/H)] with galactic mass, in conflict 
with the other studies reported here. 

The Graves et al. (2007) results are based on a large, ho- 
mogeneous data set and take advantage of new advances in 
stellar population modelling. Furthermore, Graves & Schi- 
avon (2008) have tested their method on integrated spectra 
of MW globular clusters. For a small number of clusters, 
[C/Fe] and [N/Fe] have been measured in spectra of indi- 
vidual cluster stars. The latter authors demonstrate that 
the values of [(C/Fe)] and [(N/Fe)] measured in integrated 
cluster spectra using their method match the average abun- 
dance determinations from individual cluster stars to within 
±0.04 dex. We therefore consider their abundance measure- 
ments to be reliable. However, the reader should bear in 
mind that there is not yet a consensus on C and N abun- 
dances in elliptical galaxies. 

2.2 Uncertainties in deriving [C/Fe] and [N/Fe] 
from integrated spectra 

The key results of this analysis consist of a comparison of the 
observed values of [(C/Fe)] and [(N/Fe)] from Graves et al. 
(2007) versus those obtain from chemical evolution models. 
In this context, it is important to review the uncertainties 
in the derivation of [(C/Fe)] and [(N/Fe)] using the method 
of Graves & Schiavon (2008) . We list these in order of their 
likely importance, from most to least important. 

i) Modelling [{N/Fe) ] using a CN molecular absorption 
feature. The abundance modelling process used in Graves 
et al. (2007) uses a CN molecular absorption feature to de- 
termine [(N/Fe)], after [(C/Fe)] has been fixed based on a 
C2 feature (see Graves & Schiavon 2008 for details). This 
means that the measured value of [(N/Fe)] depends on the 
accuracy of the [(C/Fe)] determination and therefore has a 
larger associated uncertainty than the other abundance ratio 
measurements. In Graves & Schiavon (2008), this method is 
tested in globular clusters and is shown to reproduce the 
average [(N/Fe)] determined from individual cluster stars to 
within 0.04 dex. This suggests that the modelling process 



is accurate, however the test has only been performed on 
two globular clusters to date. Recently, Toloba et al. (2009) 
suggested that the N abundance does not strongly correlate 
with mass by using a NH feature in the near UV, which may 
be a cleaner discriminant of N abundance. However, Toloba 
et al. did not attempt any conversion between the index and 
the actual N abundance, therefore we cannot use their data 
as a further constraint for our models. 

ii) Uncertainty in the [O/Fe] abundance. The [(C/Fe)] and 
[(N/Fe)] stellar abundance ratios are measured using the 
C24668 and CNi molecular absorption features. Both of lines 
are sensitive to the abundance of O, which competes with 
the formation of C2 and CN. Thus higher O abundances 
will inhibit the formation of these molecules, resulting in 
weaker absorption features for fixed values of [C/Fe] and 
[N/Fe]. Unfortunately, O is currently unmeasurable in spec- 
tra of unresolved stellar populations. The abundances pre- 
sented here are computed assuming solar O/Fe abundance 
ratio (i.e., [O/Fe] = 0.0). However, early type galaxies are 
known to be enhanced in a-elements and likely have super- 
solar [O/Fe]. Accounting for the likely super-solar [O/Fc] 
will have the effect of increasing the measured stellar val- 
ues of [(C/Fe)] and [(N/Fe)], since larger C and N abun- 
dances will be needed to produce the same absorption line 
strengths. As will be shown, the measured values of [C/Fe] 
and [N/Fe] are relatively high already. Also, recent evidence 
from MW bulge stars suggests that [O/Fe] may be ~solar 
for metal-rich bulge stars, despite their super-solar values of 
[Mg/Fe] (Fulbright et al. 2007). 

iii) The effects of rotation on evolutionary tracks. There is 
some indication that stars with fast rotation have bluer evo- 
lutionary tracks than non-rotating stars of the same mass 
(Meynet, private communication), which will affect the stel- 
lar population analysis of these galaxies. However, while the 
temperature differences should have a substantial impact 
on age measurements, they should have little effect on the 
abundance ratios. This can be seen in the comparison of 
solar-scale and a-enhanced evolutionary tracks in Schiavon 
(2007) and Graves et al. (2007), where turn-off temperatures 
and age measurements are strongly affected but abundance 
ratios are not. 

iv) The effects of emission infill. The galaxies in Graves et 
al. (2007) have been selected to have no detectable emission 
lines in their spectra, based on the strong emission lines at 
Ha and [OII]A3727. Thus the effects of H<5 emission are un- 
likely to substantially contaminate the CN absorption index 
in these galaxies. For galaxies in Graves et al. (2007) that 
did contain significant emission (not used in this work) , vari- 
ations in the emission infill correction resulted in differences 
of < 0.01 dex in [(N/Fe)]. 

v) Lick index zero-point offsets. The galaxy data from 
Graves et al. (2007) are not zero-point shifted to the origi- 
nal Lick/IDS system. However, this should not have a sig- 
nificant effect on the abundance determinations. Both the 
SDSS data and the Schiavon (2007) stellar population mod- 
els used to interpret them are based on flux-calibrated spec- 
tra. Schiavon (2007) demonstrates that the zero-point offsets 
between Lick indices measured in flux-calibrated spectra are 
extremely small (see Figure 2 of that work) , with the excep- 
tion of the indices Fe50I5, Fe4383, H7A, and HJa, none of 
which are used the in abundance analysis of Graves et al. 
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3 THE MODEL 

3.1 The PM04 chemical evolution model 

The adopted chemical evolution model is an updated version 
of the multi-zone model of Pipino & Matteucci (2004, PM04 
hereafter) . We calculate the evolution of element abundances 
by means of the equation of chemical evolution (see e.g. 
Matteucci & Greggio 1986 for a comprehensive discussion 
of this equation). The Initial Mass Function (IMF) <f>(m) oc 
is normalized to unity in the mass interval 0.1 — 
100M©. We adopt an IMF with x = 1.35 (Salpeter, 1955). 
We adopt the following law for the star formation rate: 

ip(t) = v • M gaa (t) , (1) 

where the star formation rate ip(t) is assumed to be propor- 
tional to the gas mass via a constant v which represents the 
star formation efficiency. We assume v = vpMm, i.e., as in 
model II of PM04 where v is an increasing function of the 
galactic mass. 

The star formation history of a model galaxy is deter- 
mined by the interplay between the infall time-scale, the 
star formation efficiency and the onset of the galactic wind 
(i.e. the energetic feedback from SNe and stellar winds). We 
assume that ip = after the development of the galactic 
wind. 

A fundamental component of this model is the detailed 
calculation of supernova (SN) explosion rates. For Type la 
SNe, we assume a progenitor model made of a C-0 white 
dwarf plus a red giant (Greggio & Renzini, 1983; Matteucci 
& Greggio, 1986). The predicted Type la SN explosion rate 
is constrained to reproduce the present day observed value 
(Mannucci et al., 2008) 

Here we adopt the same formulation for the feedback as 
in Pipino et al. (2002), to which we refer the reader. In brief, 
we consider a ~ 20% mean efficiency in energy transfer from 
the SN into the interstellar medium. In addition, we define 
the onset of the galactic wind (t gw ) as the time at which 
the energy input by SNe exceeds the gas binding energy (for 
details see PM04 and Pipino et al. 2002). The wind carries 
out the residual gas from the galaxies, thus inhibiting further 
star formation. 

In PM04, we simulated the creation of the spheroid as 
due to the collapse of either a single large gas cloud or several 
smaller gas lumps. The inclusion of gas infall makes the star 
formation rate start at a lower value than in the closed box 
case, reach a maximum, and then decrease as the gas is used 
up. 

The infall term gives the rate at which primordial gas 
is accreted by the proto-galaxy. The adopted expression is: 

{ -37 jinfall = A^in/auOe , \A) 

where Gi(t) is the mass density of element i at time t and 
Xi t i n f a ii describes the chemical composition of the accreted 
gas, which is assumed to be primordial. C is a constant 
obtained by integrating the infall law over time and requiring 
that ~ 90% of the initial gas has been accreted at t gw (in 
fact, we halt the infall of the gas at the onset of the galactic 
wind). Finally, r is the infall time-scale. 

In order to compare our predicted abundances with the 
observed ones, we must compute the mean stellar abun- 



dance of the element X ((X/H) = (Zx)), defined as (Pagel 
& Patchett 1975): 

(Zx) = ^ J ° Zx(S)dS, (3) 

where So is the total mass of stars ever born contributing 
to the light at z=0. For massive ellipticals, results obtained 
by averaging over the stellar mass are very close to those 
obtained by averaging over the stellar luminosity at z=0 (the 
difference for the models presented in this paper is typically 
less that 0.05 dex), since the indices are weighted by V-band 
luminosity (see e.g. Matteucci et al., 1998). 

3.2 General overview of the models 

We run models for elliptical galaxies with masses 10 10 Mq 
and 10 12 Mq. The basic features common to all models are 
listed in Table 1, where the input luminous mass, effective 
radius, star formation efficiency, infall time-scale and time 
of the galactic wind onset are listed in columns 1-5, respec- 
tively. Mi um is the nominal mass of the object, i.e. the mass 
of the initial gas cloud, with the infall law normalized such 
that 90% of the mass is accreted between t = and t ~ t gw . 
The mass in stars at z = is ~ 0.2 — 0.4 Mi um for all models 
and the velocity dispersion a is evaluated from the relation 
M = 4.65 • 10 s (a/km s" 1 ) 2 (i? e// /kpc) M (Burstein et al., 
1997). The effective radius R e ff is the final one, achieved 
when the collapse is over. 

Parameters such as the star formation efficiency (u), 
the infall time-scale (r), the IMF, and the fraction of binary 
systems that give rise to a SNIa explosion have been taken 
to be the same as in the best model of PM04. The latter 
authors show that these choices lead to good agreement with 
a large set of optical observables. In the present work, we 
can study for the first time the chemical evolution in much 
finer detail. The modifications to the standard PM04 nucle- 
osynthesis presented here are modest adjustments affecting 
Ca, C, and N abundances. They do not alter the satisfactory 
agreement between the observations discussed in PM04 and 
their best model predictions. 

3.3 Stellar Yields 

Although in the present work we will discuss only Mg, Fe, 
Ca, N and C, our code follows in detail the evolution of 
21 chemical elements, for which we need to adopt detailed 
stellar nucleosynthesis prescriptions. 

We first define a fiducial model (Model PM04) with the 
same nucleosynthesis prescriptions as in PM04, namely: 

(i) For single low and intermediate mass stars (0.8 < 
M/Mq < 8) we make use of the yields of van den Hoek & 
Groenewegen (1997, vdHG) as a function of metallicity, with 
the mass loss parameter tjagb dependent on the metallicity 
as follows: tjagb = 1 for Z = 0.001, tjagb = 2 for Z = 0.004, 
and tjagb = 4 for Z = 0.008, 0.02 and 0.04. A lower value 
of tjagb implies a larger yield of carbon because lower mass 
loss rates lead to longer stellar lifetimes and hence to more 
thermal pulses. As a consequence, more C is dredged up to 
the stellar surface (see Chiappini et al. 2003b). 

(ii) For massive stars (M > 8Mq) we adopt the yields 
of Thielemann et al. (1996, TNH) computed for the solar 



Table 1. Summary of model properties: physical quantities 





R eff 


V 


T 


tg W 


(M ) 


(kpc) 


(Gyr- 1 ) 


(Gyr) 


(Gyr) 


10 10 


1 


3 


o.s 


1.30 


10 12 


10 


22 


0.2 


0.44 



Table 2. Summary of model properties: stellar yields 



Model 


yields for C,N 

low- and intermediate 

mass stars 


massive stars 


PM04 


vdHG 


TNH 


Model I 


vdHG 


H07 for Z < 10~ 8 
MM02 otherwise 


Model IZ 


vdHG 


H07 


Model II 


MM02 


H07 for Z < 10 -8 
MM02 otherwise 


Model IIZ 


MM02 


H07 for Z < 0.004 
MM02 otherwise 



chemical composition. With this choice of stellar yields, the 
N from massive stars has a secondary origirjf], whereas a 
fraction of the N coming from low- and intermediate-mass 
stars has primary origin (vdHG). 

(iii) Finally, we use the yields of Nomoto et al. (1997, 
model W7) for SNIa. These are assumed to originate from 
C-0 white dwarfs in binary systems that have accreted ma- 
terial from a companion (the secondary), reached the Chan- 
drasekar mass, and exploded via C-deflagration. 

In the next sections we will also modify the base model 
to include other nucleosynthesis prescriptions. For C, N, and 
O the above fiducial model will be compared with mod- 
els computed using the stellar yields of the Geneva group, 
which include the effects of both mass loss and rotation. The 
adopted yields are the same as the ones adopted by Chiap- 
pini et al. (2006a), where details can be found. Briefly, the 
yields of Meynet & Maeder (2002, hereafter MM02) were 
adopted for metallicities larger than Z = 10~ 5 over the 
whole mass range, assuming a rotational velocity of 300 km 
s _1 . Below this value we adopted the yields of Hirschi (2007, 
H07) for massive stars, while keeping the MM02 calculations 
for low and intermediate mass stars. For SNIa we kept the 
same prescription as in PM04 model. 

The different cases studied here, with different stellar 
evolution prescriptions, are summarized below and in Table 
2. 

(i) Model I: as PM04, but the CNO yields for massive 
stars (m > 8M Q ) are as in Chiappini et al. (2006a). 

(ii) Model IZ: as Model I, except that for CNO we adopt 
the yields for massive stars at Z = 10 -8 (H07) as if they 
were valid for the entire metallicity range. The effect is to 

2 A secondary element is a chemical species which is created from 
the metals originally present in the composition of a given star. In 
the Simple Model for chemical evolution (e.g. Mattcucci 2001) it 
has been shown that the mass abundance in the gas of a secondary 
element evolves as the square of the total metallicity. 
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Table 3. Mass- weighted stellar abundance ratios 



M Lum 


[(Mg/Fc)] 


[<Fe/H>] 


[(Ca/Fe)] 


[<C/Fe>] 


[<N/Fe>] 


(M ) 












PM04 


10 10 


0.15 


0.03 


-0.07 


-0.55 


-0.19 


10 12 


0.30 


0.16 


0.04 


-0.5 


-0.18 


I 


10 10 


0.15 


0.03 


-0.07 


-0.02 


-0.05 


10 12 


0.30 


0.16 


0.04 


0.14 


-0.01 


IZ 


10 10 


0.15 


0.03 


-0.07 


0.22 


0.02 


10 12 


0.30 


0.16 


0.04 


0.34 


0.09 


II 


10 10 


0.15 


0.03 


-0.07 


-0.002 


-0.34 


10 12 


0.30 


0.16 


0.04 


0.154 


-0.19 


IIZ 


10 10 


0.15 


0.03 


-0.07 


0.09 


-0.27 


10 12 


0.30 


0.16 


0.04 


0.236 


-0.14 



increase the output yields of C and N over those of Model I. 
Although this is clearly not physically justified, we use this 
model to illustrate the effect of boosting C and N. 

(iii) Model II: in this case, we adopted the nucleosynthesis 
prescriptions for CNO used in Chiappini et al. (2006a) for 
the full stellar mass range. 

(iv) Model IIZ: as Model II, but here the H07 yields for 
massive stars are adopted for Z < 0.004. This could apply 
to systems with high star formation rates such as ellipticals, 
where the number of fast rotators could have been larger, 
extending to larger metallicities than in the MW (see De- 
cressin et al. 2007 for the same suggestion in the case of 
globular clusters). 



4 RESULTS 

In Table 3 we present the mass-weighted average stellar 
abundance ratios (cf. Eq. [3} predicted by each of the models 
described in the previous section. In the following sections, 
we analyse each of the models in detail. We first focus on 
Mg and Ca, whose abundances are well reproduced by the 
PM04 model without any further modification of the stellar 
yields. We then discuss C and N. 

4.1 Mg and Ca 

4.1.1 [(Mg/Fe)] 

The analysis of the mass-[(Mg/Fe)] relation has been com- 
prehensively addressed in PM04, Pipino & Matteucci (2006, 
2008), and Pipino et al. (2009), where a detailed compar- 
ison with the observations has been made. The Graves et 
al (2007) observational result for this abundance ratio is re- 
markably similar to previous work (e.g. Nelan et al. 2005, 
Smith et al. 2007, Bernardi et al. 2003), therefore we do not 
repeat the analysis here. We just recall that an overabun- 
dance of Mg relative to Fe is the key indicator that galaxy 
formation occurred before a substantial number of Type la 
SNe could explode and contribute to lower the [Mg/Fe] ratio 
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(for the time-delay model, see Matteucci 2001). In addition, 
the [(Mg/Fe)] ratio in the cores of ellipticals increases with 
galactic mass (Worthey et al. 1992; Weiss et al. 1995; Nelan 
et al. 2005). This relation seems to be already in place at 
redshift 0.4 (Ziegler et al 2005). 

In order to account for this trend in the star formation 
time-scale there are at least three possibilities that have been 
discussed in the literature. One involves the loss of resid- 
ual gas via galactic winds that are initiated earlier in the 
most massive objects (the inverse wind picture, see Mat- 
teucci 1994). Another possibility is to assume an initial mass 
function (IMF) which becomes systematically flatter with 
increasing galactic mass. A selective loss of metals could 
also be the cause for the increase of [(Mg/Fe)] in more mas- 
sive galaxies (see Matteucci et al. 1998). We consider the 
first possibility to be the best motivated, and we will refer 
to it as chemical downsizing. Maiolino et al. (2008) have pre- 
sented observational evidence for chemical downsizing as far 
back as z ~ 3. 

Recently Pipino et al. (2009) showed that a more phys- 
ically motivated value for v gives results in excellent agree- 
ment with the observations (and with PM04). In particular, 
following Silk (2005), they argue that a short (10 6 - 10 7 
yr) super-Eddington phase can provide the accelerated trig- 
gering of star formation needed. The subsequent quenching 
of star formation is accomplished by the SN energy input 
and results in the usual black hole mass-spheroid velocity 
dispersion relation (Magorrian et al. 1998). 

The observational scatter in the relation between 
[(Mg/Fe)] and galactic mass can be entirely explained as 
intrinsic scatter. Local effects, such as variations in the SN 
feedback efficiency of the order of a factor of two with re- 
spect to the best model case, can induce or delay the onset 
of the galactic wind and thus contribute to setting the final 
value for [(Mg/Fe)]. 

On the other hand, such a relation cannot be explained 
by a sequence of dry-mergers that create a massive spheroid 
starting from low-mass a-enhanced building blocks with- 
out violating other fundamental constraints (Pipino & Mat- 
teucci, 2008). 

The PM04 best model predicts d[(Mg/Fe)]/dloga = 
0.37, in excellent agreement with the reported value of 0.36 
by Graves et al. (2007). This abundance ratio was used by 
PM04 as a crucial constraint in order to set all the relevant 
physical quantities reported in Table 1. The same values are 
adopted here so that Models I, IZ, II and IIZ also obey to 
this fundamental constraint. 

In the following sections we study other abundance ra- 
tios and, when a discrepancy is found between model pre- 
dictions and the observations, the problem is addressed by 
comparing model predictions computed with different stellar 
yield prescriptions. 

4.1.2 [(Ca/Fe)] 

The observations show that, although Ca belongs to the 
group of the a elements, the strength of the observed 
Ca lines appear to follow [(Fe/H)] instead of [(Mg/H)] 
(Worthey, 1998; Trager et al., 1998; Saglia et al. 2003). 
In particular Thomas et al. (2003) suggested [(Ca/Mg)] 
=—0.15, and Saglia et al. (2003), after a detailed analysis 
of several possible sources of error, claimed that the Ca de- 



pletion in ellipticals is real. Such a result is confirmed by 
Graves et al. (2007), who find the [Ca/Fe] ratio to be nearly 
solar over the entire galactic mass range. 

The PM04 calculations already suggested [(Ca/Mg)] = 
—0.152 and [(Ca/Fe)] = —0.03 (these values are for the core 
of a 10 n M G PM04 galaxy) and are in excellent agreement 
with the above mentioned observations. Here we add that at 
either lower or higher galactic masses the predicted ratio is 
close to the solar value with a very mild relation with mass. 
A linear regression fit of the mass weighted stellar abundance 
ratio as a function of the stellar velocity dispersion returns 
d[(Ca/Fe)]/d log cr = 0.27 for PM04 (and hence the same 
value for Models I, IZ, II and IIZ), steeper than the slope 
of 0.13 reported by Graves et al. (2007). We notice that 
if we adopt Woosley & Weaver (1995) yields for massive 
stars we obtain d[(Ca/Fe)] jd log a = 0.18, which is in better 
agreement with the observed slope but produces [(Ca/Fe)] 
values offset about 0.15 dex higher than observations at a 
given mass. 

As shown in PM04, this result can be explained simply 
by SN yields. In fact, the a-elements exhibit different de- 
grees of enhancement with respect to Fe in the results of 
the chemical evolution model (see Francois et al. 2004). This 
is due, from a theoretical point of view, to the different de- 
gree of production of each element in Type II and la SNe. In 
particular, Ca and Si show a lower overabundance relative 
to Fe than O and Mg, since they are also produced substan- 
tially in Type la SNe. In fact, the mass of Ca ejected during 
a SNIa explosion in the model W7 (Nomoto et al., 1997) 
is ~ 0.012Mq, whereas the contribution of a typical Type 
II SNe (averaged over a Salpeter IMF in the mass range 
10 - 5OAf ) is ~ 0.0058M Q (see Table 3 of Iwamoto et al., 
1999). While some fine tuning of the stellar nucleosynthesis 
could improve even further the agreement with observations, 
we stress that PM04 explains the behaviour of the Ca abun- 
dance as a function of mass without requiring changes in the 
adopted chemical evolution model parameters. 



4.2 C and N 

Figs.Q]and [2]show our model predictions for the [(C,N/Fe)]- 
mass relations compared with the observed ones by Graves 
et al. (2007) (represented by pentagons connected by a solid 
line. The PM04 Model (triangles joined by a solid line) pre- 
dicts both [(C/Fe)] and [(N/Fe)] to be below solar and to 
show no dependence on galactic mass, at variance with the 
observations. The other models (i.e. I: dotted line; IZ: dashed 
line; II: long-dashed line; IIZ: long-dashed dot line), while 
keeping the same star formation history as PM04, show dif- 
ferent trends in abundance ratios with mass, reflecting the 
different assumptions about the nucleosynthesis. Most of the 
modified models lead to an improvement over the predictions 
of PM04. Thus, for the first time we can use local elliptical 
galaxies as a test bench for stellar nucleosynthesis studies. 
In the following sections will we discuss the C and N results 
in further detail. 



4-2.1 [(C/Fe)] in the integrated stellar spectrum 

We first focus on the predictions of our new models con- 
cerning the [(C/Fe)] abundance ratio in Fig[T] Model I (dot- 
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Figure 1. Mass- weighted [(C/Fc)] as a function of galactic veloc- 
ity dispersion predicted by PM04 (triangles linked by a solid line) 
and Models I (dotted line), IZ (short-dashed line), II (long-dashed 
line) and IIZ(short dash-dotted line) compared to the observed 
relation found by Graves et al. (2007, pentagons linked by a solid 
line). 



tial additional C enhancement from massive stars but is not 
based on any physical motivation. 

Finally, Model IIZ matches the [(C/Fe)]-mass relation 
quite well. In this case, the C contribution from fast rotators 
is extended up to Z — 0.004 and is therefore larger than in 
Model II (where fast rotators exist only up to Z = 10 -8 , see 
Section 3.3). 

Overall, it is encouraging that all of the models which 
include the effects of rotation and mass loss in massive 
stars (Models I, IZ, II, and IIZ) produce reasonable agree- 
ment with the data, showing the correct slope of the mass- 
[(C/Fe)] relation, with zero-point offsets of the order of only 
±0.1 dex from the data. There are variations depending on 
the various yield tables adopted, but these variations are 
relatively small. This suggests that modelling of the nucle- 
osynthesis of C in massive rotating stars might approach 
a consensus soon and therefore that C abundances will be 
a powerful new tool for constraining galaxy star formation 
histories. 

It should be kept in mind that the transformation of 
Lick indices into mean stellar abundances carries uncertain- 
ties, as discussed in Section 2.2. These uncertainties are more 
likely to affect the zero-point calibration of the abundances, 
while the relative abundances and therefore the abundance 
trends with a should be fairly robust. Finally we notice that 
both our model predictions and the data set of Graves et 
al. (2007) are at odds with the recent claim (Clemens et al. 
2006) of a decrease in the [(C/H)] value in stars as function 
of the galactic mass. 



ted line) predicts d[(C/Fe)]/dlog a = 0.4, in good agree- 
ment with the value of 0.5 reported by Graves et al. (2007, 
solid line) . However, this model has a zero-point of the mass- 
[(C/Fe)] relation ~0.1 dex smaller than the observed one. 

The large difference between Model I and PM04 is due 
to the larger contribution of massive stars to the carbon pro- 
duction in the case of the Geneva group yields with respect 
to the TNH ones. The larger C yields of the Geneva group 
are found not only at very low metallicities (due to rotation, 
as explained before) but also at solar metallicities (see Fig.l 
of Chiappini et al. 2003b for a comparison of the C yields 
of the two groups at Z = Zq). In fact, the mass-loss rate, 
which increases with metallicity, contribute to increase the 
overall C yields (see MM02). 

A similar result is found for Model II (long-dashed line), 
which differs from Model I only in the C yields for interme- 
diate mass stars (vdHG in Model I and MM02 in Model II). 
The difference between these two models is expected to be 
small since, as shown by Chiappini et al. (2003a) the yields 
of C from vdHG and MM02 for low- and intermediate-mass 
stars are similar. In fact, although MM02 yields do not in- 
clude the third dredge up contribution, they obtain a similar 
effect thanks to rotation and mass-loss (see Chiappini et al. 
2003a for details). 

Model IZ (dashed line) predicts: d[(C/Fe)]/dlog a = 0.3 
with a larger zero-point. In this case the model galaxies are 
more C-enhanced than the observed ones. This is not sur- 
prising, as Model IZ assumes fast rotators to exist at all 
metallicities and not only below Z=10 -8 as done in Chi- 
appini et al. (2006a). This illustrates the effect of substan- 



4-2.2 [(N/Fe)] in the integrated stellar spectrum 

Concerning N, models I and IZ (see Fig symbols as in 
Fig. [TJ, substantially increase the zero-point of the mass- 
[(N/Fe)] relation with respect to PM04, hence improving 
the agreement with the observations. However, both models 
predict a rather flat slope: d[(N/Fe)]/dloga = 0.1 (a factor 
of 5 flatter than the value indicated by the SDSS galaxies). 
We remind the reader that Models I and IZ assume primary 
N production in massive stars, whereas in PM04 the N from 
massive stars is only secondary. 

On the other hand, Model II is able to predict a slope 
(d[(N/Fc)]/dlog<7 = 0.4) for the mass- [(N/Fe)] relation, in 
much better agreement with Graves et al.'s findings, but it 
does not improve upon the original PM04 model in terms of 
the absolute value. In fact, [(N/Fe)] is under-predicted by 
0.4 dex (i.e. more than a factor of 2) in this case. The only 
differences between models I and II are the adopted yields 
for low- and intermediate-mass stars. Model I was computed 
with the yields of vdHG, which predict larger quantities of N 
than the yields of MM02 (adopted in Model II). These large 
N quantities are responsible for both the larger zero point 
and the flatter behavior of the mass- [(N/Fe)] relation. When 
assuming that the fast rotators play a role up to metallicities 
Z = 0.004 (Model IIZ), the [(N/Fe)] zero point is increased 
by ~0.1 dex, which is still not enough to account for the 
observed values. 

In the case of Models II and IIZ, the differences between 
observed and predicted values at a given galactic mass are 
probably too large to be explained in terms of either un- 
certainties in the calibration of the Lick system (see Section 
2.2), or in the actual measurements. This offset is larger than 
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Figure 2. Mass-weighted [(N/Fe)] as a function of galactic ve- 
locity dispersion predicted by PM04 (triangles linked by a solid 
line) and Models I (dotted line), IZ (dashed line), II (long-dashed 
line) and IIZ (short-dashed dot line) compared to the observed 
relation found by Graves et al. (2007, pentagons linked by a solid 
line). 

the 3cr measurement error quoted in Graves et al. (2007) for 
the observed zero-point of the [(N/Fe)]-<r relation. 

It is worth noting that, unlike in the case with [(C/Fe)], 
there is substantial disagreement in the predicted values of 
[(N/Fe)] between the four models, both in the predicted val- 
ues of d[(N/Fe)]/a!log a and in the zero-points (at the level 
of 0.4 dex). This indicates that the various nucleosynthe- 
sis yield predictions for N vary significantly between groups 
and are highly sensitive to the mass and metallicity range 
over which rotation is assumed to be relevant. The nucle- 
osynthetic predictions for N are therefore significantly less 
secure than those for C. 

Also, as discussed in Section 2.2, stellar [(N/Fe)] abun- 
dance ratio measurements in the observed galaxies are more 
uncertain than the other abundance determinations. 



4-2.3 N abundance in the gas: model ellipticals 

Since the predicted values for [(C,N/Fe)] are the result of 
the convolution of the temporal evolution of [C,N/Fe] in the 
gas with the star formation rate, it is important to study 
the gas abundances in order to understand the differences 
between the models presented in the previous section. In 
Fig. [3] we plot the evolution of [N/Fe] in the gas phase of 
the 10 10 M© and the 10 12 M© model galaxies (left and right 
panels, respectively) as a function of [Fe/H]. The dotted lines 
mark relevant times during the course of star formation in 
the galaxy, including the metallicity at the time when the 
SFR reaches its peak, the metallicities at which the SFR is 
within a factor of two from the peak value, and the onset of 



the galactic wind (labelled as t gw ). This event sets the time 
(hence the metallicity) at which the star formation stops. 
This illustrates the values of [N/Fe] at which most of the 
stars in each galaxy are formed. 

First of all, a comparison of Figs [2] and [3] shows that, 
even if N is more enhanced in more massive (hence more 
metal rich) galaxies, this does not imply that N has a sec- 
ondary origin! It is true that models in which N is mostly 
secondary (e.g., PM04) have gas phase [N/Fe] abundances 
which increase steeply with [Fe/H], as shown in Figure [3] 
(PM04 is the solid black line). However, this does not neces- 
sarily translate into a strong [(N/Fe)]-mass relation (see Fig- 
ure [2| . In contrast, models which produce a better match to 
the observed [(N/Fe)]-mass relation (e.g., Model IZ, dashed 
blue line) , include substantial primary N production in mas- 
sive stars and show high gas-phase [N /Fe] abundances at low 
[Fe/H] (Figure 0. 

By comparing the left panel of Figure [3] (10 10 Mq case) 
to the right one (1O 12 M0 case), it is clear that the more mas- 
sive galaxy features a higher [N/Fe] ratio in the gas during 
the phases when the star formation proceeds at its maximum 
rate. Stars form out of this enriched material until the onset 
of the galactic wind at time t gw . The difference between the 
1O 1O M0 galaxy and the 10 12 M© galaxy is the smallest for 
the PM04 model, hence we predict a rather flat slope in the 
[{N/Fe)]-mass relation. The difference is largest in models 
II and IIZ, which therefore exhibit the steepest slope in the 
[(N/Fe)]-mass relation (Fig. 2). 

This result is quantified in Fig.Ufor Model IIZ: the nor- 
malized distribution of stars as a function of [N/Fe] peaks at 
higher values of the abundance ratio and features a smaller 
tail at [N/Fe] < -1 in the 10 12 M Q model galaxy (full his- 
togram) than in the 10 10 Mq galaxy (empty histogram). 
Note that the values of [(N/Fe)] derived by means of Eq.(3) 
and shown in Table 3 correspond to the mean value of dis- 
tributions like the one shown in Fig. U The mean does 
not always coincide with the peak value (i.e. the mode) of 
the distribution. Given the asymmetry of the distributions, 
the mean, in fact, tends to be ~0.1-0.2 dex lower than the 
mode. By comparing the distribution of stars as a function 
of [N/Fe] for different models, we can also explain the dif- 
ferences in the predicted [(N/Fe)] at a given mass (Section 
4.2.2). For instance, model IZ features a distribution which 
peaks at higher [(N/Fe)] and has a larger mean value than 
model PM04. Finally, we note that all the models shown in 
Fig. 3 (with the exception of models IZ and IIZ) exhibit a 
[N/Fe] ratio below solar during most of the galaxy evolution. 

Intermediate mass stars (i.e., M < 8Mq) do not con- 
tribute metals to the ISM until ~ 30 Myr after the beginning 
of star formation. This corresponds to Z = 5 • 10~ 4 in the 
10 10 Af o galaxy and Z = 3 ■ 10~ 3 in the 10 12 M Q galaxy, 
where the evolution is faster. Due to the very short time- 
scale of galaxy formation (see Table 1), the lowest stellar 
mass which can contribute to the chemical enrichment of the 
model galaxies is ~ I.QMq in the 1O 1O M0 case and ~ 2Mq 
in the 1O 12 M case. The differences in the N production at 
low metallicities (i.e. [Fe/H] and [O/H] < —2) will have 
little effect in the final [(N/Fe)] ratio because very few stars 
are formed at these metallicities. 
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Figure 3. [N/Fe] as a function of [Fe/H] in the gas for a 10 10 Mq (left panel) and a 10 12 Mq (right panel) elliptical galaxy. Different 
lines represent the different models in Table 2: PM04 (solid line) and Models 1 (dotted line), IZ (dashed line), II (long-dashed line) and 
IIZ (short-dashed dot line). Vertical dotted lines indicate the metallicity where the SFR reaches its peak, as well as when it rises to and 
drops by a factor of two from this value. Also shown is the metallicity at the onset of the galactic wind (labelled as t gw ), indicating the 
time (hence the metallicity) at which the star formation stops. This illustrates the values of [N/Fe] at which most of the stars in each 
galaxy are formed. 




Figure 4. Normalized number of stars as a function of [N/Fe] for model IIZ in the case of the 10 10 Mq model galaxy (empty histogram) 
and of the 10 12 Mq model (filled histogram). 



4-2.4 N abundance in the gas: high redshift progenitors 

Ideally, one would use gas phase abundance ratios as shown 
in Fig [3] to further refine the model. Unfortunately, no such 
measurements exists in low-redshift ellipticals, as they are 
almost devoid of gas. In order to perform such an exer- 
cise we must resort to high-redshift objects, such as the 
Lyman Break Galaxy (LBG) MS1512-cB58 (Pettini et al. 
2001, Teplitz et al. 2000). This is currently the brightest 



LBG known, because it is gravitationally lensed. The LBG 
is at z — 2.7672 and has a luminous mass of ~ 10 10 M© 
with a star formation rate of ~ 40Moyr -1 (Pettini et al. 
2001). We use this galaxy because, from the analysis of its 
properties and its chemical abundance pattern, Matteucci & 
Pipino (2002) suggest MS 1512-cB58 to be a low-mass young 
elliptical experiencing its main episode of star formation and 
galactic wind. 
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Figure[5]shows the evolution of the [N/O] abundance ra- 
tio as function of [O /H] in the gas for our model galaxy with 
Mi, — 1O 1O M0, for each of the nucleosynthesis cases studied 
in the previous sections. In this figure, the point with error- 
bars represents the LBC{f]. The models shown in Fig. [5] have 
star formation rates similar to the one observed in MS1512 
cB-58. The solid line represents the original PM04 model. 
The [N/O] abundance ratio rises quite steeply with metal- 
licit)|j because, in the PM04 model, the N from massive stars 
is assumed to be a secondary element. The predictions of the 
PM04 model (which is an improved version of the chemical 
evolution model presented by Matteucci & Pipino, 2002) are 
within la of the observed value. 

Model II (long dashed line) performs equally well, 
exhibiting an early (i.e. [O/H] < —3) plateau with val- 
ues around [N/O] = —1, followed by a decrease due to 
the metallicity-dependent yields of massive stars. Later, 
intermediate-mass stars start contributing to the chemical 
enrichment (once [O/H] > —1 dex in the case shown in 
Fig. [5} . Such a contribution starts to be important around 
[O/H] = —0.5, where a new rise in the [N/O] ratio is ob- 
served. 

The detailed behaviour of the curve predicted by Model 
II is still uncertain, primarily due to the following reasons: 

(a) There are still no complete grids of stellar yields for 
models with rotation. Here we are combining the calcula- 
tions of MM02 for 300 km s _1 with the more recent calcula- 
tions of H07, which includes the fast rotators (500-800 km 
s _1 ). Grids with intermediate values of rotational velocities 
at different metallicities are currently not available. More- 
over, there are small differences between MM02 and H07 
calculations, such as different amounts of convective over- 
shooting, that prevent a completely self-consistent approach. 

(b) The contribution to N by Super-AGB star^f] could be 
important (Siess 2007). This is not included here as no stel- 
lar yield calculations for these objects are available at the 
moment. These stars would play a role at earlier times than 
low- and intermediate-mass stars and could contribute to 
fill up the valley seen in Fig. [5] (also contributing to increase 
the median value of [(N/Fe)] in Fig. E}. (c) In MM02, the 
intermediate-mass stars were not followed up to the end of 
their evolution and the contribution of hot bottom burning 
could be missing (but see Chiappini et al. 2003b). 

Despite of the above caveats, the important point here is 
that Model II predicts a high N/O ratio at early phases that 
later decreases and then increases again. The amount of the 
decrease in N/O at intermediate metallicities is still uncer- 
tain but it clearly shows that some metallicity-dependence of 
the N yields in massive stars is necessary in order to bring 
[N/O] back up to levels consistent with the Lyman Break 
galaxy. 

On the other hand, Models IZ and IIZ, which are built 
to enhance C and N in order to fit the stellar ratios in local 

3 Note that the N abundance is inferred from absorption lines 
(Pettini et al., 2001), whereas the O abundance comes from emis- 
sion line measurements (Teplitz et al., 2000) 

4 O contributes 60% of the total metallicity in this system. 

5 Super-AGB stars are massive enough to ignite C in a degenerate 
core but do not proceed into neon burning and beyond. They are 
thus an intermediate case between intermediate-mass and massive 
stars. 



ellipticals, predict a large N abundance with respect to O 
at low [O/H]. This is because we let massive stars produce 
N with the same stellar yields as if they were in the low- 
est metallicity regime. The sudden downturn of Model IIZ 
at [O/H] ~ —0.5 is due to the adopted threshold metallicity 
for including the fast rotators. This drop is much less impor- 
tant than what happens in Model II, where the fast rotators 
contribute during a much shorter period of the evolution 
of these galaxies. In other words, in Model IIZ the drop is 
just shifted to larger metallicities because we allowed fast 
rotators to contribute for longer times (up to larger metal- 
licities) . 

Model I (dotted line) is an intermediate case between 
PM04 and Model II. At variance with the fiducial PM04 
model, it features an early enhancement of N owing to the 
production from low-metallicity massive stars, as in Model 
II. It also shows a steeper rise around [O/H] = —0.5 due to 
the fact that in this case, vdHG nucleosynthesis for inter- 
mediate mass stars is producing more N (the hot bottom 
burning contribution). 

While Model I and IIZ are still marginally consis- 
tent with the data from Pettini et al. (2001, triangle in 
Fig-EJ, Model IZ — which better matches the [(N/Fe)]-mass 
relation — predicts a [N/O] ratio more than a factor of three 
larger than the observed LBG value at [O/H] = —0.35. 
However, the gas-phase N abundance inferred from emis- 
sion lines in the same galaxy is [N/O] ~ —0.5 (Teplitz et 
al. 2000, full square in Fig. [5}, more than three times higher 
than the [N/O] inferred from the absorption line measure- 
ments of Pettini et al. (2001). Pettini et al. argue that the N 
abundance inferred from emission lines may be substantially 
overestimated. Nonetheless, we can interpret the factor of 
three difference in the N abundance between emission and 
absorption lines as an estimate of the uncertainty in the gas- 
phase measurements. Using the Teplitz et al. (2000) emission 
line abundance brings Model IZ into excellent agreement 
with both stellar- and gas-phase abundances (full square in 
Fig- 0). Unfortunately, they do not provide errorbars, there- 
fore we cannot quantitavely estimate the level of agreement 
with our model prediction. 

It is important to stress that a comparison with a single 
high-redshift object, which my be a young elliptical caught 
in the act of formation, may not be appropriate. More data 
on gas phase abundances in high redshift objects are needed 
to confirm the tension emerging with respect to abundances 
measured in stars and to reconcile emission- and absorp- 
tion- line measurements of the gas-phase abundance. How- 
ever, along with the observations of DLA (which also show 
systematically lower N/O ratios than the stars in the MW, 
again suggesting a Z-dependency of the N yields in massive 
stars), such a comparison demonstrates the existence of a 
dichotomy between stellar abundances (which tend to favor 
higher N abundances) and gas phase abundances. 

Our analysis demonstrates a potential inconsistency be- 
tween the stellar N abundances observed in local ellipticals 
and the gas-phase abundances observed in a likely progeni- 
tor of these galaxies at higher redshift. The stellar [(N/Fe)] 
abundance ratio appears to require substantial N enrich- 
ment due to both massive and intermediate-mass stars (e.g., 
Model IZ), but this model produces an [N/O] abundance ra- 
tio that is too high to match the LBG data from Pettini et 
al. (2001). This discrepancy disappears if the [N/O] abun- 
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Figure 5. [N/O] as a function of [O/H] in the gas for the different 
10 10 M Q models: PM04 (solid line) and Models I (dotted line), 
IZ (dashed line), II (long-dashed line) and IIZ (short-dashed dot 
line). The triangle with error-bars is the observed value for the 
Lyman Break Galaxy MS1512-cB58 (Pettini et al 2001), whereas 
the solid square is the estimate from Teplitz et al. (2000). 



dance ratio of MS1512-cB58 is in fact [N/O] = -0.5, as 
determined by Teplitz et al. (2000). 

Another possible solution to this discrepancy is that 
the [(N/Fe)] ratios inferred from stellar spectra are over- 
estimated by about 0.3 dex (this is the amount necessary 
to bring it into agreement with model prediction and the 
Pettini et al. LBG data at the same time). The various un- 
certainties in the [(N/Fe)] measurements were discussed in 
Section 2.2. An offset of 0.3 dex is larger than can easily 
be accounted for in this way, but is certainly not impos- 
sible, given the uncertainties in the modelling process, the 
paucity of good test data, and the fact that these constitute 
the first attempts to quantify [(N/Fe)] in unresolved stellar 
populations. 



5 CONCLUSIONS 

In this paper we have analysed the mass- [Ca,C, N/Fe] ratio 
relations observed in local spheroids by means of a chemi- 
cal evolution model which successfully reproduces the mass- 
[(Mg/Fe)] and the mass-metallicity relations (the PM04 
model) . Our theoretical predictions are compared to the ob- 
served relations inferred from integrated spectra of early- 
type galaxies in SDSS by Graves et al. (2007). We find that: 

• The [(Ca/Fe)]-mass relation is naturally explained by 
the PM04 best model without any further assumption 
needed. In particular, the under-abundance of Ca with re- 
spect to Mg can be attributed to the relative contributions 



of SNe Type la and Type II to the nucleosynthesis of these 
two elements. 

• The standard nucleosynthesis prescriptions adopted by 
PM04 (i.e. where rotation and mass loss are not taken into 
account) lead to chemical evolution models in which C traces 
Fe and N behaves as a secondary element (i.e., there is no 
primary nitrogen contribution from massive stars). Instead, 
the abundance ratios measured for the stellar population in 
local spheroids suggest that both C and N behave more like 
a-elements (e.g., Mg) in the sense that the [(C,N/Fe)] abun- 
dance ratios are super-solar and increase with the galactic 
mass. 

• We show that in order to reproduce the [(C/Fe)]- and 
[(N/Fe)]-mass observed relations it is necessary to consider 
new nucleosynthesis prescriptions, which take into account 
both the effects of rotation and mass-loss (Meynet et al. 
2009). These two important physical processes have already 
proven to be crucial for our understanding of the abundance 
ratios in metal-poor stars in the Galactic halo (Spite et al. 
2005, 2006). Indeed, the high N/O, C/O and low 12 C/ 13 C 
ratios observed in stars with [Fe/H] < —2.5 suggest the ex- 
istence of fast rotating stars in the very early Universe (Chi- 
appini et al. 2006a,b, 2008). 

• The various different prescriptions for the effects of ro- 
tation and mass loss on the nucleosynthesis of C considered 
here all produce results that are consistent with the data and 
with each other. This indicates that the nucleosynthesis of 
C is now understood well enough and hence C abundances 
can be used to place constraints on chemical evolution in 
galaxies. The same level of consensus is not found for N 
abundances from the various models. This suggests that the 
nucleosynthesis of N is more sensitive to the differences in 
stellar evolution models and that the existing yield tables 
for N are perhaps less reliable than those for C. 

• Once these new stellar yields are implemented in the 
PM04 chemical evolution model, a remarkable agreement 
with the observations is attained (within the observational 
uncertainties) both for the slope and the zero-point of the 
mass-[(C/Fe)] relation (Model II). This shows that signifi- 
cant amounts of C must be produced in massive stars (and 
not only in intermediate mass stars) in order to reproduce 
the [(C/Fe)] abundances of massive early type galaxies. 

• The same model reproduces the steep slope observed for 
the mass- [(N/Fe)] relation in the local universe. This model 
also provides predictions for N/O ratios in the gas phase 
abundance of young spheroids. The models considered here 
cannot simultaneously predict the high values of [(N/Fe)] 
observed in local stellar populations and the low values of 
[N/O] observed in the gas-phase of LBG MS1512-cB58. A 
factor of ~ 0.3 dex is required to resolve this difference. 
There are significant uncertainties in both the stellar and 
gas-phase measurements of N, as well as the N yields, all of 
which may contribute to the observed discrepancy. 

In summary the PM04 model, updated to include new 
nucleosynthetic yields from stellar models which account 
for stellar rotation and mass loss, can account for the 
[(Mg,Ca,C/Fe)]-mass relations observed in local spheroids. 
Discrepancies may remain in the case of N. However, our 
best model (Model II) might be completely reconciled with 
observations if emission-line gas phase N abundance of MS 
1512-cB58 is considered, rather than the N abundance de- 
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termined from absorption-line measurements. Another way 
to reconcile the discrepancy is if the [(N/Fe)] in the galaxies, 
inferred from Lick indices, is overestimated by ~ 0.3dex. If 
this turns out to be the case, our results suggest that fast 
stellar rotation at low Z and large mass-loss rates at higher 
metallicities in massive stars have an important impact on 
the chemical evolution of early-type galaxies, as has already 
been shown for the Milky Way. On the other hand, if the 
[(N/Fe)] ratios inferred for the SDSS galaxies are correct 
and if the absorption-line gas phase N abundance of MS 
1512-cB58 are typical of young spheroids, a tension between 
observed gas and stellar abundance ratios remains, suggest- 
ing that further processes not envisaged in the present work 
should be taken into account. However, the uncertainties as- 
sociated with both stellar and gas-phase N abundance mea- 
surements, as well as the uncertainties in the nucleosynthesis 
of N do not allow us to draw firm conclusions. 

Finally, we note that, thanks to large sample sizes, high 
quality data and improved stellar population models, we are 
able for the first time to use elliptical galaxies as a test 
bench for stellar nucleosynthesis studies. At the same time 
we present several predictions for the evolution of the N /Fe 
and N/O abundance in the gas that might be tested on high- 
redshift proto-galaxies. 
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